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Abstract 

The consistency condition is tested within the particle-particle random-phase approximation 
(RPA), renormalized RPA (RRPA) and the self-consistent RPA (SCRPA) making use of the 
Richardson model of pairing. The two-particle separation energy is calculated in two ways, namely 
as the energy of the first addition mode, which adds two particles to a core with N particles, and 
as the energy of the first removal mode, which removes two particles from the N + 2 particle system 
to get back to the same iV-particle core. The corresponding transitions generated by the pairing 
operators are also calculated. It is found that the results obtained in these two ways of calcula- 
tions are close to each other only at large values of particle number N and/or small interaction 
strength. At N < 10 for a given value of the interaction strength, the discrepancy between the 
results obtained in two ways of calculations within the SCRPA is much smaller than those given 
by the RPA and RRPA. 
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I. INTRODUCTION 



The random-phase approximation (RPA) has been widely used in the theoretical study of 
nuclei within the valley of /5-stability. The success of the RPA is mainly based on the use of 
the quasiboson approximation (QBA), which considers fermion pairs as boson operators, just 
neglecting the Pauli principle between them. As a result, a set of linear equations, called the 
RPA equations, was derived, which reveals the physics of collective excitations generated by 
the RPA boson-like modes. The simplicity of the RPA equations allows a feasible treatment 
of a number of many-body problems, which would be computationally intractable otherwise. 
However, this approach suffers a drawback: It breaks down at a certain critical value of the 
interaction's parameter, where the RPA yields imaginary solutions. The reason of this well- 
known RPA instability is the violation of Pauli principle within the QBA. 

In /3-stable medium and heavy nuclei, the QBA is a good approximation, and the RPA is a 
very powerful tool for the description of several important quantities such as the ground-state 
and excited-state energies, electro-magnetic transition probabilities and their distribution, 
transition densities, etc. The first-order diagram expansion beyond the mean field treated 
within the RPA includes significant effects of two-body correlations beyond the mean-field. 
However, with reducing the particle number, the concept of collective excitations, which are 
described by the RPA modes, becomes less and less firm. The ground-state correlations 
(GSC) which are left beyond the RPA become stronger in light systems. This feature makes 
the validity of the QBA, and therefore of the RPA itself, questionable in the systems with 
small particle numbers. 

Several approaches were developed to take into account the GSC beyond RPA in a simple 
way such as to restore the Pauli principle among the fermion pairs, from which the RPA 
operators are constructed. The popular one, known as the renormalized RPA (RRPA) Q, 
l2L Li \M- includes the expectation value over the ground state of the diagonal elements of 
the commutator between two fermion-pair operators, neglected in the QBA. In this way the 
Pauli principle is approximately taken care of. The inclusion of GSC beyond RPA within 
the RRPA eventually renormalizes the interaction in such a way that the collapse of RPA 
is avoided, and the RRPA has solutions at any value of the interaction's parameter. This 
scheme has been tested in several exactly solvable schematic models such as the Lipkin 
model for ph RRPA Q, |^, the SO(8) model for the renormalized quasiparticle RPA 
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and the Richardson model for the pp RRPA and self-consistent RPA (SCRPA) 

UU- The 

results of these studies showed that the difference between the energies of the first excited 
state given by the pp RPA, pp RRPA, and SCRPA increases noticeably with increasing the 
interaction parameter and/or decreasing the particle number. The SCRPA solutions are the 
closest ones to the exact solutions. These results raise a serious question of the applicability 
of the RPA to the calculations for realistic light-neutron rich nuclei. As a matter of fact, a 
close agreement between the RPA results and the experimental data for light neutron-rich 
nuclei can be considered really reasonable if the corrections due to GSC beyond RPA are 
negligible, which justify the validity of the QBA 8]. In the opposite case such agreement 
may be even illusory. 

This concern is the motivation of the present work, whose goal is to test the consistency 
of the pp RPA and its renormalizations, the pp RRPA and SCRPA, at different particle 
numbers as the interaction strength varies. The reason of choosing the pp RPA comes 
from one of its major merits, which allows the straightforward calculation of an important 
quantity in the study of unstable nuclei, namely the two-particle separation energy. Indeed, 
the pp RPA uses the addition and removal phonon operators to create the ground states 
and excited states of the systems with N + 2 and N — 2 particles, respectively, from the 
ground state of the iV-particle system. Therefore the two-particle separation energy of the 
(N + 2)-particle system can be calculated either as the energy of the lowest excited state 
generated by the addition mode, which adds two particles to the core with N particles, or as 
that generated by the removal mode, which removes two particles from the core with N + 2 
particles. A consistent theory should give identical results in either way of calculation. By 
analyzing the results obtained in two ways of calculations within all three approaches, namely 
the RPA, RRPA, and SCRPA in their application to an exactly solvable schematic model, 
namely the Richardson model Q], the present work may shed light on the consistency of 
these approaches, in particular, in their application to systems with small particle numbers 
(light systems). 

The paper is organized as follows. The brief outline of the pp RPA, RRPA, and SCRPA 
applied to the Richardson model is presented in Sec. |Hj As the present paper deals only 
with pp RPA and its renormalization, the prefix pp will be omitted hereafter. The results of 
numerical calculations are analyzed in Sec. IIHI The paper is summarized in the last section, 
where conclusions are drawn. 
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II. RPA, RRPA AND SCRPA WITHIN THE RICHARDSON MODEL 



A. The Richardson model 

The Richardson model, considered in the present paper, consists of Q doubly-fold equidis- 
tant levels, which interact via a pairing force with parameter G. The model Hamiltonian is 
given as ja ll | 

n n 
H = J> - X)N, -GJ2 P}Pi , (1) 

where the particle-number operator iVj and pairing operators P-, Pi are defined as 

N t = c\c t + cUc, , P} = 4^, Pi = {P}) ] - (2) 
These operators fulfill the following exact commutation relations 

[P u P}] = 8 ij (\-N i ) , (3) 

[N h P}] = 25^], [Ni, Pj] = -MyPj . (4) 

The single-particle energies of the equidistant levels are defined as e« = ie with i running 
over all Q levels. The present paper considers only the p/i-symmetric case. This means that, 
in the absence of interaction (G =0), the lowest Q/2 levels are filled with N = Q particles 
(two particles on each level). Numerating particle (p) and hole (h) levels from the levels 
nearest to the Fermi level, the particle and hole energies are equal to e p = e(0/2 + p) and 
€h = e(0/2 — h+ 1), respectively, with p (h) = 1, . . . , Q/2. The Fermi level ep is defined as 

e F = [e(n + l)-G]/2 . (5) 

There are several methods of solving exactly the pairing problem described by Hamiltonian 
namely the Richardson method jjj, the infinite-dimensional algebras in Ref. Q], and 



the direct diagonalization of Hamiltonian ((TJ in the Fock space 11]. 



B. RPA 



The outline of the RPA is given in detail in Ref. It uses the addition and removal 
Pl\ phonon operators to create the ground states and excited states of the systems with N+2 



and N — 2 particles, respectively, from the ground state of the iV-particle system. Applied 
to the Richardson model, these addition and removal phonon operators are given as 



and 



respectively, where 



4 = E^M-E^ 



Ql = Pi , Q h = -Pi 



(6) 



(7) 



(8) 



The RPA assumes the validity of the QBA, which replaces the expectation value of the 
commutation relations [Qi, Qj] in the RPA ground state |iV, 0) of iV-particle system with 
that obtained in the Hartree-Fock (HF) ground state |HF), namely 

(N,0\[Qi,Q]]W,0) ~ (HF|[Qj, Q]]|HF) = 6 tj , {i,j) = (p,p') ,(h,ti) , (9) 

since (HF|A^|HF) = 0, and (HF|JV h |HF) = 2. Using the QBA ©, one can easily verify that 
the expectation values of the commutation relations for addition and removal operators in 
the RPA ground state are of bosonic type, i.e. 

(N,0\[A fM ,Al,]\N,0) = 6^ , (N,0\[R x ,R{,]\N,0)=6 xy , (10) 

if the X and Y amplitudes obey the following normalization (orthogonality) conditions 



while the closure relations 



(ii) 



£4 A ^-£x^ = o 

A fi 

guarantee the following inverse transformation of Eqs. © and ((7j) 



Ql 



Qi 



E4 A) ^+E^ (M) 4 



(12) 



(13) 



The RPA equations are obtained by linearizing the equation of motion with Hamiltonian 
(JIJ and operators A* and W. The matrix form of the RPA equation for the addition mode 
is 

A b\ ( XM 



-B C 



E,. 



( iM 

y(/0 



1 



where the submatrices A, B, and C are found as 



A^ A =(HF|[Q„, [J/,Q»,]]|HF) = 2 





G" 


H) 


+ 2 



(14) 



(15) 



R RPA 



^hh' 



-<HF|[Q P , [^,Q h ]]|HF) = G 
-<HF|[Q h , [ff,Q*,]]|HF) = -2 



<W + G 



(16) 



(17) 



C. RRPA 



The QBA Q neglects GSC beyond RPA. The latter make the expectation value of N p 
deviate from zero and that of Nh deviate from 2. In order to take this effect into account, 
the RRPA considers the following renormalized addition and removal phonon operators 

4 = E ~ E - ^d ftt = J- X^Ql - Y P (X) Q P , (18) 

p h hp 

respectively , with the abbreviation 
O] 



— o 



(A) = (A,0|A|A,0) , i = p,h 



(19) 



where the correlated ground state \N,0) is defined as the vacuum of the operators and 
Tlx, i.e. 

a^\n7o) = n x \NTo) = o . (20) 



The GSC factors (A) are the expectation values of the operators 

D p = l-N p , D h = N h -l, 



(21) 



in the correlated ground state |iV, 0) ([20)1 . They are introduced in the definition (|18|) to 
preserve the ground-state expectation value of the commutation relation Q. Indeed, using 



Eq. © and definition (JHJ), one obtains the exact commutation relations for operators Qi 
and Q\ in the form 

[Qpi Q p '] — bpp'Dp , [Qh, Qh>] = Shh'Dh • (22) 

Using this exact relation (|22|) and the definition (|2U|) . one can see that the renormalized 
addition and removal operators satisfy the boson commutation relations in the correlated 
ground state (j2"Uj) 

([^,4']> = V> ([n x ,n{,]) = 5 xy , (23) 

if the amplitudes X and Y satisfy the same RPA orthogonality conditions while the 
same RPA closure relations (|12j) guarantee the following inverse transformation of Eq. (|18|) 



Ql = v w E x p" )a I + E y P (A) ^ 



E4 A) ^ + E^4 



(24) 



instead of Eq. (fTB^l . 

The RRPA equations are obtained in the same way as that for the derivation of the RPA 
equations, assuming the factorization 



(25) 



as well as neglecting all the expectation values (Q p ,Q p ), (Q P Qh), and (Q h Qh')- The RRPA 



submatrices obtained in this way for the addition modes have the form 

i rrpa _ 2 



A itrtl 



1\ G 
r \"~2 + 2 



5 PP > - GJ (D p )(D p/ ) 



(26) 



R RRPA 

tiph 



GJ(D p )(D h ) 



(27) 



^RRPA 



<W + Gy/{D h ){D h ,) . 



(28) 



1\ G 
2j + 2 

They are different from the RPA submatrices (JT3J) - (JT7J) by the square roots (Dj) (Dj), 
which renormalize the interaction. In the absence of GSC beyond RPA the expectation 
values (Di) become 1, and the RPA submatrices are recovered. 
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The equations for the expectation values (_Dj) are obtained following Refs. in the 
form 

w = tt^?- {d ^tt^v' (29) 

which are exact in the present model because of the exact relation 

N t = 2P}P i . (30) 



For the details see Ref. [6|. The set of Eqs. 1(1 ^1 . pfi jl - (|23 Jl . and forms the closed 
set of RRPA equations, which are non-linear due to the presence of the backward-going 
amplitudes Y in the expressions for the GSC factors [Dj) (|29|) . In numerical calculations 
this set is solved self-consistently by iteration. Knowing (-Dj), one obtains the occupation 
numbers fi from Eq. ({21} as 

f P = \{l-{D p )), f h =1(1 + <£>*)) , (31) 

so that within the RPA (Dj = 1) the HF occupation numbers f p = and fh = 1 are 
recovered. 



D. SCRPA 



The only difference between the RRPA and SCRPA is that the latter includes the so- 
called screening factors, which are the expectations values (Q p iQ p ), (Q P Qh), and (Q h Qh>) 
neglected within the RRPA. They are derived using the inverse transformation ()24|) and the 
definition of the correlated ground state (|2U|l. and have the form ((JO]: 

A 

= (OlOj) = -<PiP P > = -(P}Pk) = ^(D p )(D h )J2Xh X) Y p {X) , (33) 

A 

(QhQh>) = V / WWE y /^' ) = < P i^> - *wW , (34) 

with (p+p ft ) = vw(^)E4 A) 4' ] • 

A 
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The SCRPA submatrices are then given as 



A pp , = ([Q p , [H,Ql,]}) 



1\ G 



e p 



+ 



G 



h" 



V ^4== . (35) 



^ = -<[Q p , [ir,g fc ]]> = G 



_[DjM_ 



C hh , = -([Q h ,[H,Ql\\) 



(36) 



G 



G 



/i" 



<W + G } DhDv) =. (37) 



In principle no assumption on the factorization (J25)) is required within the SCRPA. However 
the factorization (|25j) does simplify greatly the calculations, whose results turn out to be 
quite close to those obtained without assuming (|23jl [6]. The equations for the GSC factors 
are the same as Eq. (J2HJ). 



E. Two-particle separation energy and pair transition 

1. Consistency in terms of two-particle separation energy 

Unlike the particle-hole RPA, where the solutions with negative eigenvalues have no 
physical meaning, each set of the (particle-particle) RPA, RRPA, and SCRPA equations 
has m + n solutions, all of which correspond to physical states. From them m (= N/2 in 
the present case) solutions correspond to the system with + 2 particles, whose excitation 
energies related to the ground-state (g.s.) energy of the N-particle system are E^ = 
£( N + 2 ) _ gW) w ith /i = 1, . . . , m. The other n (= N/2) solutions have the eigenvalues equal 
to — E x = — [^g^? — with A = 1, . . . , n. They correspond to the eigenstates of the 

system with N — 2 particles with the eigenvalues E\. Therefore the first eigenvalues E p= \ of 
the (A^+2)-particle system and E\ = \ of the (N — 2)-particle one are the ground-state energies 
of these systems related to that of the A^-particle core. This means that they correspond to 
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the two-particle separation energies S£ ) (N + 2) in the (N + 2)-particle system, and S^ p '(N) 
in the core, which are defined as 

(N + 2) = E, =1 = - £g> , S£ (N) = E X=1 = £g} - £^ r 2) • (38) 

As such, in an exact or consistent theory, it follows from Eq. f)38jl that 

S^(N + 2) = si;\N + 2) . (39) 

This equation is called the consistency condition hereafter, which means that the two ways 
of calculating the two-particle separation energy as the energy of the addition mode, which 
adds two particles to the TV-particle core, and that of the removal mode, which removes two 
particles from the (iV + 2)-particle system to get back to the same core, should give the same 
result. The expressions (|58j) were used in Ref. ^| to calculate the two-neutron separation 
energies for N c + 2 and N c nuclei, namely 12 ' 14 C, 10 ' 12 Be, and 9,11 Li, whose neutron cores 
consist of N c = 8 particles. 



2. Pair-transition probability 



The pair transition is generated by the operator 

k=J2 (p! + Pi) = + Qp) - + Qh) ■ 

i=p,h p h 



(40) 



By using the inverse transformation (|24jl . one obtains the square | (N, 0\Aij,K\N, 0) | 2 of the 
matrix element of the operator K, which corresponds to the transition probability from the 
ground state of the iV-particle system to the one-phonon state of the (N + 2)-particle 
system within the RRPA and/or SCRPA, in the following form 

B^N -^N + 2) = KivTol^livTo)! 2 = T, \[W)X^ - V vWn W 2 • (41) 



In the same way, the transition probability from the ground state of the iV-particle system to 
the one-phonon state | A) of the (N — 2)-particle system is proportional to | (N, 0\TZ\K\N, 0) | 2 , 
which is given as 

2 



Bx (N- 2 <- N) = | {N, 0\K x K\N,0) p 



(42) 
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The corresponding expressions within the RPA are recovered by setting (Di) = 1 in Eqs. 
(j4*Tj) and (|32j) . From the consistency condition discussed in the preceding section it follows 
that 



which means that the probabilities of the pair transition between the ground states of N- 
and (N + 2)-particle systems should be the same in either direction. 

In the limit with e p> \ = e p= \ — e/2 and eh>\ = £/i=i = — e/2, the Richardson model is 
reduced to the well-known two-level model, in which each of the two levels is iV-fold and the 
Fermi lev el (151) i s located at e-p = —G/2. This degenerate model has been studied in detail 



in Refs. mm, Il5j |. The RRPA solutions of this model yields the following two-particle 
separation energies [See, e.g, Eqs. (73) and (74) of Ref. |l5| ] 



Sg\N + 2) = y/(e + G)[e-G{(D)N-l)] , Sg{N + 2) = y/(e + G)[e - G((D)N + 1)] . 



The RPA result is recovered by setting (D) = 1. Result P4*j) analytically shows that, for this 
degenerate two-level model, the consistency condition (jSUjl is asymptotically fulfilled only 
when N > 1. 

III. ANALYSIS OF NUMERICAL CALCULATIONS 

The calculations were carried out within the Richardson model with the level distance 
e = 1 MeV at several values of the particle number N as the pairing parameter G varies. 

Shown in Fig. ^are the two-particle separation energies S$(N + 2), which are obtained 
within the RPA, RRPA, and SCRPA according to Eq. (|3%j) at N = 2. 10, 50, and 100. 
For a comparison, the exact results, available for N = 2 and 10, are also shown. For iV = 
50 and 100 the size of the matrix makes the exact diagonalization infeasible so that the 
exact solutions are not available. In Ref. it has been analyzed in detail that the RRPA 
and SCRPA extend the solutions far beyond the point where the RPA collapses, and that 
the SCRPA predictions are the closest ones to the exact solutions. Therefore the details of 
these features are not repeated here. This figure shows that, for iV < 10, the consistency 
condition (}3T?j) is clearly violated in all the approximations under consideration. The violation 
is particular strong within the RPA at small N. It is seen that (N + 2) > (N + 2) 



B^ X {N ^N + 2) = B X=1 (N <- N + 2) , 



(43) 




(44) 
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FIG. 1: Two-particle separation energies S^n (N+2) at several values of N as functions of pairing 
parameter G. The dashed, dotted, and dash-dotted lines represent results obtained within RPA, 
RRPA, and SCRPA, respectively, for which the thin lines denote S-C\n + 2), while the thick lines 
stand for S^(N + 2). The exact results for N = 2 and 10 are shown by the thick solid lines. 

within the RPA and RRPA, but S$(N + 2) < S$(N + 2) within the SCRPA. At a 
given value of G, the SCRPA predicts the smallest difference between (N + 2) and 
S2 n \N + 2), while the RPA gives the largest one. One can see how (N + 2) gets closer 
to S$(N+2) as N increases so that Eq. is perfectly restored at N = 100. The small 
difference between (N + 2) and (N + 2) and the better agreement with the exact 
result, obtained within the SCRPA even at small N, show that the SCRPA is the most 
consistent approximation among three approximations under consideration. At the same 
time, the quite large difference between (N + 2) and (N + 2) obtained within the 
RPA at small shows that the RPA is, in principle, an inconsistent theory when applied 
to light systems. This observation means that, GSC are important and cannot be neglected 
in calculations for light systems. 



12 




FIG. 2: Occupation numbers fh of five hole levels nearest to the Fermi level (h = 1,. . . , 5) for the 
systems with N = 10 (a) and 100 (b). The open circles, squares, and full circles denote the RPA, 
RRPA, and SCRPA results, respectively. The triangles in (a) are the exact results. The dashed, 
dotted, dash-dotted, and solid lines connecting these discrete points are drawn to guide the eye. 
The results in (a) and (b) are obtained by using G = 0.33 MeV and 0.17 MeV, respectively. 

In order to see the origin of the effect due to GSC, the occupation numbers fh of five hole 
levels, which are located nearest to the Fermi level, are plotted in Fig. |2]for A^ = 10 and 100. 
They are obtained by using the values of G equal to 0.33 MeV and 0.17 MeV for A^ = 10 
and 100, respectively. These values are just below the collapsing points of the RPA, which 
are G crit = 0.34 MeV and 0.18 MeV for N = 10 and 100, respectively. The results in Fig. 
121 (a) shows that the QBA, on which the RPA is based, is no longer a good approximation 
when A" is small. Indeed, the QBA means that fh = 1, while the value of fh, obtained for 
the level just below the Fermi level, is around 0.7 within the RPA. At large N, e.g. A" = 
100, the QBA is much better fulfilled as the occupation number fh=\ is around 0.85, while 
those for all the other 49 hole levels are larger than 0.95 [Fig. |21 (b)]. This figure also shows 
that the effect of GCS, once included within the RRPA and SCRPA, makes the RRPA and 
SCRPA phonon operators much closer to the ideal bosons since the values of fh predicted by 
the RRPA and SCRPA for the hole level just below the Fermi level are around 0.9 and 0.95, 
respectively, for A^ = 10, and larger than these values at larger N. The comparison of the 
occupation numbers predicted by the RPA, RRPA, and SCRPA again shows that the QBA 
is best satisfied within the SCRPA. It is interesting to notice that, while the SCRPA predicts 
a two-particle separation energy much closer to the exact value for A^ = 10 as compared to 



13 



10 



8 



-h 6 
4 

CQ 



(a) N = 2 



12 



1 1 1 r jr^r 

(c) N = 50 / // 

« // 

/// 
/# 

// 



0.2 0.4 0.6 0.1 
16 



(b) N = 1 



12 

CQ 



i i 
i i 
i 

•• i 
i 
i 

I ! 



1 o 

n 20 r 

15 
10 



0.05 0.1 0.15 0.2 0.25 



(d) N = 1 00 



I ; I 

i // 

/ // 

I 1 1 

I ! I 

I I i 

UJ 



/// 

/// 



0.1 0.2 0.3 0.4 0.05 0.1 0.15 0.2 

G (MeV) G (MeV) 



FIG. 3: Squares of matrix elements corresponding to the pair-transition probabilities as functions 
of pairing parameter G at several values of N . The notations for dashed, dotted, and dash-dotted 
lines are as in Fig. ^ The thick and thin lines denote B fl= i(N — > N + 2) and B\ = \(N <— N + 2), 
respectively. 

the RPA [Fig. Q(b)], the agreement between the SCRPA prediction and exact values for 
the occupation number jh=\ is not as good as that offered by the RRPA. 

The pair-transition probabilities, shown in Fig. |3J increase with increasing G as one might 
expect. In a similar way as that for the two-particle separation energies, the condition (J4*H|) 
is strongly violated within the RPA for small N. Combining the results on Figs. ^ and |31 
it is seen that for N < 10, neither condition nor ()4Hj) holds. Like the restoration of the 
consistency condition (J3nj) for two-particle separation energies, as N increases the condition 
(jlHJ) for pair-transition probabilities is also gradually restored to become well-fulfilled at 
N = 100. 
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IV. CONCLUSIONS 



The present paper tested the consistency of the RPA, RRPA, and SCRPA within the 
Richardson model. The consistency condition under consideration here is the requirement 
that the two-particle separation energy s!£'(N + 2) obtained as the energy of the first 
addition mode, which adds 2 particles to the iV-particle system, should be the same as 
S2 n \N + 2) obtained as the energy of the first removal mode, which removes 2 particles 
from the N + 2-particle system. The results of calculations show that this consistency 
condition is well fulfilled only for sufficiently large values of the particle number iV (> 50 
in the present model), i.e. in medium-mass and heavy systems. For light systems (N < 
10) it is found that, among the three approximations under consideration, the RPA, which 
has the largest deviation from the exact results at large G, strongly violates the consistency 
condition. The SCRPA agrees best with the exact solutions and also is the most consistent 
approximation in the sense that the above-mentioned consistency condition is well satisfied 
already at iV > 10 even at G above the RPA collapsing point. 

The results obtained in this test indicate that GSC beyond RPA becomes quite important 
in light systems, especially in the region close to the RPA collapsing point. They invalidate 
the QBA, which is assumed for the pair operators within the RPA. Therefore, in principle, the 
RPA predictions of the two-neutron separation energies in light neutron-rich nuclei cannot 
be reliable without taking the effect of GSC beyond RPA into account. From the three 
approximations considered in the present paper, the most consistent one, which includes this 
effect for light systems, turns out to be the SCRPA. Although the present test was conducted 
within a schematic model with doubly-fold equidistant levels, it is unlikely that this general 
conclusion on the inconsistency of the RPA would change when applied to realistic light 
nuclei. As a matter of fact, the recent calculations of the two-neutron separation energy 
5*2n in 12 ' 14 Be by using the Gogny interaction have shown that the RRPA shifts the result 
by around 14% for 12 Be, and the consistency condition (jHOj) is not fulfilled (N = 10) for 
both the RPA and RRPA 16] . An easy way to see if the RPA needs to be renormalized is 
to estimate the occupation numbers (f3Tj) using the backward-going amplitudes Y obtained 
within the RPA. If there is a strong deviation of f p (fh) from (1), this means that the 
validity of the RPA is deteriorated. In the above-mentioned calculations for Be-isotopes e.g., 
the occupation number of the last occupied shell amounts to around 0.8 within the RPA, 
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and increases to around 0.95 within the RRPA, showing the necessity of using the SCRPA, 
or at least the RRPA, instead of the RPA. 
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